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Abstract 

In this letter, the Kondo magnetic effect is studied for the XX Z spin chain where 
the impurities are coupled to the edges of the system. The Hamiltonian of our model 
can be constructed from the transfer matrix. It is exactly solvable and the exchange 
constants between the bulk and the boundary impurities are arbitrary. The finite 
size corrections to the ground state energy are obtained for the magnetic impurities 
and the boundary condition. The specific heat and the susceptibility contributed 
by the impurities are derived and the Kondo temperature is given explicitly by the 
use of the Luttinger- Fermi liquid picture and the Bethe ansatz method. 
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It is well known that the spin dynamics of the Kondo problem is equivalent to the dynamics of the spin chain with 
the magnetic impurities Jl]] . Much progress has been made recently on the impurity models based on the methods of 
renormalization-group techniques conformal field theory [^-|| and the integrability investigations. The quantum 
inverse scattering method and the Bethe ansatz technique have been the powerful tools to study the integrable 
impurities within the framework of the quantum spin chains. 

Magnetic impurities play an important role in the strongly correlated electron systems, especially for one dimensional 
systems. And the impurity usually destroys the integrability when it is embedded into the 'pure' quantum chain. 
Then it is a challenging problem to deal with the impurity effects in the strongly correlated electron systems for the 
quantum integrable models without losing the integrability. Now several integrable impurity models have been found 
|6|-|9| besides the magnetic impurities in the noninteracting metals such as the Anderson impurity model and the 
multichannel Kondo problem p|,^0[-p2|. 

As a quantum mechanical model of magnetism, the Heisenberg spin chain turned out to be very fruitful and the 
integrable impurity problem for this model was first considered in Ref.[13] and generalized to arbitrary spin in Ref. 
[14]. Besides the natural understanding of completely integrable quantum systems, the quantum method of the inverse 
problem provides us also a very useful technique to construct the impurity model which preserves integrability. Indeed, 
this scheme has been adapted in Ref. [6-8]. But, for the periodic boundary systems, some unphysical terms must 
be present in the Hamiltonian to maintain the integrability, though they may be irrelevant when the impurities are 
embedded in the system. Notice that the fact that the impurities cut the one-dimensional system into the 'piece' 
when they are introduced and (then) the open boundary systems are formed with the impurities at the ends of the 
systems. So the integrable impurity models |l5|-|l8|] can be studied with the use of open boundary conditions. In recent 
years, various integrable models on the open chains with boundary terms have been investigated (See, for example, 
references [19-22]). In Ref. [23] the open Heisenberg chain with impurities has been discussed and the integrability is 
given in Ref. [24] when the two impurities are coupled to the XX Z chain. In the present letter, we will discuss the 
effects of the magnetic impurities in the XXZ spin chain and the finite size correction contributed by the impurities. 
Our starting point is the following Hamiltonian 

JV-l 

+j(^ + « + ^{«c) (i) 

where the constants J, c and 77 satisfy the relation J = sinh 2 rj cosh c/ (sinh 2 rj — sinh 2 c) . The Pauli operators a^' y ' z 
act non-trivially only in the jth quantum vector space and d^ ,y ' z , d x h v,z are the spins of magnetic impurities with 
s = 1/2 . This system can been constructed from the transfer matrix [p5| 

U (A) = Tr T [K+ (A) T T (A) K~ (A) f T (A)} (2) 

where 

T T (A) = L Tb (A + c) L tN (A) • • • L T2 (A) L rX (A) , 
f T (A) = L Tl (A) L T2 (A) • ■ ■ L tN (A) L Tb (A - c) . 

The L operator is given by 

4 

L rn (A) = ^2wj<T J T <g) a J n 
3 = 1 

with 

W\ = W2 = — sinh?7, 
W4 — W3 = sinh A, 
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W4 + Ws — sinh (A + rf) 



for the subindex n = 1, 2, • • • , N , a, b being the quantum vector spaces and the impurity sites at both of the ends. 
The symbol Tr denotes the trace in the auxiliary space t in C 2 . The <7 J are the Pauli operators for j — 1, 2, 3 with 
<jj being the identity operator. The boundary reflection matrices take their forms as 

K+ (A) - I, K- (A) = L Ta (A + c) L Ta (A - c) 

where c is an arbitrary constant. Then the Hamiltonian (|l]) can be obtained from the logarithmic derivative of the 
transfer matrix U (A) when the spectral parameter A — > . By the use of the standard Bethe ansatz method, the 
eigenvalue of the Hamiltonian is 

sinh 77 

E = 4 sinh -q > j ■ 



r 



cosh (2Xj) — cosh 77 



2J \ 

TV — 1 H cosh 77 

cosh c / 

for the Hamiltonian (jlj) with the Bethe ansatz equation 

cosh (Aj - §) s inh 2 (Aj + c + §) 



cosh (Aj + 2) sinh 2 (Aj + c - 

sinh 2 (Aj - c + 2) J s inh (Aj + f) 
'sinh 2 (Aj- C -§) jsinh(Aj -§) 

M 



y-r sinh (Aj- - Afc + 7?) sinh (Aj + A fc + 77) 
^ sinh (Aj — Afc — 77) sinh (Aj + Afc — 77) 



where j = 1, 2, • • • , M with M being the number of down spins. Let pspTf l 



$ (A, ^ ) = 2 tan" 1 



7 

cot — tanh A 



1 - F 2l -f 



By taking the logarithm of the Bethe ansatz equation (||), we have that 

2$ (Aj + c, |) + 2$ (Aj - c, I 



+ (2JV+1)$(A J( 2) -*(A is r) 

47r/j + £ {$ (Aj - A fc , 7 ) + $ (Aj + Afc, 7 )} (4) 



fc=i 

where 77 = ij . Ij is an integer or the half-odd integer. Under the thermodynamic limits we can write down the Bethe 
ansatz equation as the form 

f-OC 

dfia (fi) $' (A - /j, 7) + 2na (A) 
, 1 N 



where 
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cr N (A) = 2$'(a + c,^) +2$'( 



A — c, — 



-*'(A,r) + #' (a, J 



The a (A) is the distributed function. The a N (A) is the effect of the open boundary condition and the magnetic 
impurities, ft will give us the finite size corrections introduced by the two impurities coupled to both of the ends of 
the system. Fourier transformation on the equation (|J) gives us the distributed function as the form 

o-(A) 



2 7 cosh^ 2N(2ttY 

z2\u sinh ( 7rw ) ° N M 



due 

2 sinh (ttuj — -fuj) cosh yyw) 



with 



^ N . . 47rsinh (7ro; — 70;) 
cr (u) = 



sinh (7tcj) 

sinh (7raj — 2Tu>) 1 



f + 4 cos (2cw) 



sinh (jru) — ju>) J 



where the second term with the factor 1/ (2N) denotes the finite size correction due to open boundary impurities. 
Therefore, the distributed function has the form 



f f 
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CT(A) ^27^A + ^ a/(A) 



with the finite correction: 



°7 ( A ) = tt / dwe 



1 r 8^(^-2^) 



27r J-oo sinh (7ro; — 7^) cosh (70;) 

1 f 



7r ( A + c ) ^ p^l, tt(A-c) 

f 



7 cosh n ^~ CJ 7 cosh 



2 7 cosh 2^' 

' 7 



ft is owing to the magnetic impurities coupled to the ends of the spin chain. The ground state energy per site for our 
system is 

e . 2 r°° d\ 

2 sin 7 



N J -00 cosh (ttX) [cos 7 — cosh (2 7 A)] 

f / 2J 
N \ cosh c 



1 / 2J 
+ cos 7 + — — f 



+2 / dX- ff /( A ) sin2 T 



cos 7 — cosh (2A) 

It is obvious that the ground state energy corresponding to the bulk is exactly same as the periodic case ^8|. Its 
finite correction is, based on the boundary condition, 

E - I sin I f°° duj sinh [(^ - 2r)u] 1 
J_ 00 sinh (ttuj) cosh (70;) 

dx- sin27 



—00 



cosh (ttX) [cos 7 — cosh (27 A)] ' 
The contribution given by the magnetic impurities has the following form: 
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Ei = 



dX 



4 sin 2 7 



cosh (tt\) [cos 7 — cosh (27A — 2c)] 



2 sin 7 



sinh 2 c 



- sin 2 7 



Obviously, when c — 0, our model degenerates to the XX Z spin chain with the free boundary condition. 

In the following section, we shall study further the physical properties of the impurity model using the Luttinger- 
Fermi theory. We can regard 2irIj/N in equation (^) as the momentum of the jth quasi-particle. Similarly as in the 
periodic boundary case, the energy E takes its minimized value when A — > 00 . Then the Fermi velocity of the system 
can be defined as 



Vf 



£ 'd W 



2trt (A) 



where the symbol Ed denotes the dressed energy of the model (See reference [29]). Now we set Of = cr* + with the 
relations 



<Ti (A) 



7 cosh 



7T(A + C) 

7 



7 cosh 



tt(A-c) ' 



a b (A) =-± 
Z7T 



,-i2Au 



oc 

1 



sinh (7tcj — 2ro;) 
sinh (ttw — 70;) cosh (70;) 



27 cosh 



7TA 



describing the finite corrections of the distributed function contributed by the magnetic impurities and open boundary 
condition, respectively. The state density of the system at the Fermi surface can be expressed as 30 - 32 13] 



AT (00) 



1 



1 



1 <7j(A) 

2irv F [' ' 2iVcr (A) 
1 oib(A)' 



2A^ a (A) 



A— >oo 



Therefore, the low temperature specific heat contributed by the magnetic impurities has the expression: 

= 1 7CQsh^ ^ 
3N sin 7 

where T denotes the temperature of the system. The correction to the susceptibility is 

Xi = Icosh^^xo 

where xo is the susceptibility in the bulk. The Kondo temperature for our system is given by 

N 1 



T k = 



2tt cosh 2£ 

7 



since the Kondo temperature is nothing but the effective Fermi temperature in the local Fermi-liquid theory [p3[ . 
Notice that the Kondo temperature ~ iVcos -1 (ttc/j) when c — > ic. This is similar to the result for the XXX 
Heisenberg spin chain and supports also the result obtained earlier in Ref. [34], that there is a crossover from an 
exponential law to an algebraic law when the coupling constant between the host spins and the impurities changes 
for the Kondo temperature. 

In the above discussion, the impurity contributions are obtained in the thermodynamic limits since the Bethe 
ansatz equations can be written down as the integral equations which may be solved by a Fourier transformation. It 
is an interesting problem to discuss the impurity effects for the finite lattice case. In fact, in Ref. [35] de Vega and 
Woynarovich have given a method to calculate the leading-order finite size corrections to the ground state energy of 
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the model which is soluble by the Bethe ansatz. The situations for the XXZ Heisenberg chain and the Hubbard 
chain are also studied []36| , [37||2g| l by using the finite size scaling technique |3^|4C||. The impurity effects in different 
sectors of the present model and its critical properties are also the interesting subjects for further discussions. 

To conclude, we have studied the solvable magnetic impurity model within the framework of open boundary XXZ 
Heisenberg chain where the impurities are coupled to the ends of the system. The Hamiltonian of the system with 
impurities can been derived from the transfer matrix. The finite size correction to the ground state energy is obtained 
for the contribution of the magnetic impurities. With the help of the Luttinger-Fermi description we get the specific 
heat and the susceptibility due to the impurities. The Kondo temperature is also given explicitly. 
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